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AIR FORCES AND MOMENTS ON TRIANGULAR AND RELATED WINGS WITH SUBSONIC 
LEADING EDGES OSCILLATING IN SUPERSONIC POTENTIAL FLOW ' 


Ву CHARLES E. Жаткіхв and JULIAN Н. BERMAN 


SUMMARY 


This analysis treats the air forces and moments in supersonic 
potential flow on oscillating triangular wings and а series of 
xweptback and arrow wings with subsonic leading edges and 
supersonic trailing edges. For the wings undergoing sinusoidal 
torsional oscillations simultaneously with vertical translations, 
the linearized velocity potential is derived in the form of a power 
series in terms of a frequency parameter. This method can ђе 
useful for treatment of similar problems for other plan forms 
and for wings undergoing other sinusoidal motions. For tri- 
angular wings, as many ferms of such а series expansion as 
тау be desired can be determined; however, the terms after the 
first few become cery cumbersome. 

Closed expressions that include the reduced frequency to the 
fifth power, an order which is sufficient for a large class of 
practical applications, are given for the velocity potential and 
for the components of chordwise sectton force and moment 
coeficiente. 

These wings are found to exhibit the possibility of undamped 
torsional oscillations for certain ranges of Mach number and 
locations of the axis of rotation. The ranges of these parameters 
are delineated for triangular wings. 


INTRODUCTION 


This report is concerned with the derivation of expressions 
for the velocity potential and associated forces and moments 
for oscillating triangular wings in supersonic flow. The 
boundary-value problem for the linearized velocity potential 
for an apex-forward triangular wing oscillating іп а super- 
sonic main stream may be classified, according to reference 
1. as "purely supersonic” if the leading edges of the triangle 
are outside the Масћ cone emanating from the apex of the 
triangle or “mixed supersonic” if the leading edges are inside 
this Mach cone. 

In the purely supersonic case the principle of independence 
holds; that is, the fow on the upper surface of the wing is 
independent of the flow on the lower surface and vice versa. 
Garrick and Rubinow (reference 1) have shown that the 
boundary-value problem for the velocity potential in the 
purely supersonic case can be satisfied by simple distributions 
of sources with local strength proportional to the local pre- 
scribed normal velocity of the wing. Nelson (reference 2) 
has treated the oscillating triangular wing for this case to the 
third power of the frequency. 


In the mixed supersonic case the principle of independence 
does not hold. Boundary-value problems for lift-producing 
wings in this case can be satisfied by distributions of doublets; 
the relation between doublet strength and normal velocity 
of the wing is, however, in general, not sumple. The deter- 
mination of this relation requires the solution of an integral 
equation that employs the potential of а time-dependent unit 
doublet as kernel and limits of integration that depend on 
Mach number and wing plan form. 

For treatment of problems that involve boundary condi- 
tions that are independent of time—such as constant angle 
of attack, constant rate of pitching, and so forth—the doublet 
potential, employed as kernel of the integral equation, is 
considered independent of time and in these cases the integral 
equations for triangular wings can be solved by a straight- 
forward process. 

For treatment of problems of oscillating wings, however, 
it is necessary to employ, as the kernel of the integral equa- 
tions, а doublet potential that varies harmonically with time 
&nd in this case the solution of the integral equation, gener- 
ally, becomes very cumbersome. If the doublet potential or 
kernel is expanded in terms of the frequency of oscillation, 
however, use can be made of knowledge of solutions of inte- 
gral equations for problems that are independent of time to 
obtain &n expanded form of solution for & wing undergoing 
harmonic oscillations. Such & procedure was demonstrated 
in treatments of rectangular wings in references 3 апа 4. 

Тће purpose of the present report is to make use of the 
expanded form of the velocity potential to obtain the forces 
and moments, based on the first few terms of this potential, 
for a rigid triangular wing performing vertical and pitching 
sinusoidal oscillations in mixed supersonic fow. Although 
as many terms of the expanded potential as may be desired 
can be obtained after the first few terms, the process becomes 
very cumbersome. The flow normal to the leading edge is 
subsonic but the flow normal to the trailing edge is considered 
to be supersonic. This latter consideration implies that the 
potential derived for triangular plan forms may be used to 
calculate the aerodynamic forces and moments for other plan 
forms that may be formed with the triangular wing by cut- 
ting the trailing edges so that they lie ahead of the Mach 
cones emanating from their foremost points. | | 

Other approaches to the solution of the problem of oscillat- 
ing triangular wings have been given by Robinson (reference 
5), Haskind and Falkovich (reference 6), and by Stewartson | 


1 Supersedes МАСА TN 2457, “Air Forcesand Momentson Triangular and Related Wings With Subsonie Leading Edges Oscillating in Supersonic Potential Flow" by Charles E. 


Watkins, 1951. 
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(reference 7). Іп both references 5 and 6 formal solutions 
to the problem were obtained in terms of special systems of 
curvilinear coordinates. Robinson's solution was given in 
terms of a double summation of trilinear combinations of 
Bessel functions of the first kind with Lamé functions of the 
first and second kinds. Similarly, the solution of Haskind 
and Falkovich was given in terms of summations of Bessel 
functions of the first kind combined with elliptic integrals 
of the first and second kinds. In both references 5 and 6 the 
potentials were not reduced to useful forms for calculating 
forces and moments. 

In reference 7 Stewartson makes an interesting though 
specialized use of the Laplace transformation to develop a 
method whereby terms of the velocity potential for triangular 
plan forms, expanded as herein, can be obtained. Stewartson 
gives formulas that, except for errors presumably in printing, 
can be used to develop the potential to the second power of 
the frequency, but he omits many details in his derivation. 


SYMBOLS 
ф disturbance-velocity potentia! 
2,0,2 rectangular coordinates attached to wing 


moving in negative z-direction 





Ey rectangular coordinates used to repre- 
sent space location of doublets in 
zy-plane 

Zm function defining mean ordinates of any 
chordwise section of wing such as 
y=y as shown in figure 1 

w(x Yı, t) vertical velocity at surface of wing along 
chordwise section at Y= 

Xo abscissa of axis of rotation of wing as 
shown in figure 1 

t time 

h vertical displacement of axis of rotation 

ho amplitude of vertical displacement of 
axis of rotation, positive downward 

a angle of attack 

Og amplitude of angular displacement 
about axis of rotation, positive leading 
edge up 

ва time derivatives of А and o, respectively 

V velocity of main stream 

е velocity of sound 

M free-stream Mach number (V/c) 

B—JM*—1 | 

ധ frequency of oscillations 

- Mw 

(с == T 2 

k reduced frequency (bw/V) 

є half apex angle 

C—tan є 

Anm represents functions of w, z, and M 

„Ку functions used to denote doublet distri- 


bution functions 


constantes associated with D, depending 
on 8C 
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ДР ДЕ; constants depending on 8C 

е: constants depending on ВС and М 

2b root chord of wing 

8 semispan of wing 

c slope of ray passing through vertex of 
wing 

47,7 | dummy variables 

p density 

Ар local pressure difference 

P x section force (total force at any spanwise 
station) 

[27770277 components of section force coefficients 

Ma section moment (total moment about 
മാമ at any spanwise station) 

21,24,,24,,21; components of section moment coeffi- 
cients 

M, total component of damping-moment 
coefficient 

po ВС 

F',E' complete elliptic integrals of the first 
and second kinds, respectiv ely, with 
moduli А1 — pj 

ANALYSIS 


BOUNDARY-VALUE PROBLEM FOR THE VELOCITY POTENTIAL 


Referred fo a rectangular coordinate system moving 
forward at a uniform supersonic speed in the negative 
z-direction (see fig. 1), the differential equation for the 
propagation of small disturbances that must be satisfied by 
the velocity potential is 


1/0, рд оф п о% 
as tV Ae Orr tz 


The main governing boundary condition to be satisfied by 
the velocity potential is that the flow be tangent to the 
surface of the wing or 


(1) 


офу _ 92, Zn 
(9). —w(z, y o V em om 


(2) 
where Zn is the vertical displacement of any point of the wing. 
For the particular case of & wing independently performing 
small sinusoidal pitching oscillations of amplitude о; about 
some spanwise axis 70 and small sinusoidal vertical trans- 
lations of amplitude fo, the quantity Z, in equation (2) is 


Lm =e [ох — 9) - ho] = oft — 2o) +h (3) 


(See fig. 1 (b) for sketch showing instantaneous displacement 
of section у=.) For convenience, the frequency of oscil- 
lation of both pitehing and translatory motion is denoted by 
w. Considering these motions to occur at separate frequen- 
cies would add no difficulties to the derivation. 

Substituting the expression for Z4 (equation (3)) into 
equation (2) gives 


w(z,y,D = V a4d-à(z—29) +h (4) 
Equation (4) implies that the velocity potential may be 
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expressed as the sum of geparate effects due to position and 
motion of the wing associated with individual terms of this 
equation, namely 


ф= фа + bat bi (5) 
DERIVATION OF Ф 


In order to obtain the analytical expression for the 
potential ¢, 15 is necessary to derive only one of the subsidiary 
potentials appearing in equation (5), say фа. The other 
subsidiary potentials ¢, and ф can then be obtained from 
the derived expression for ф, by simple comparison. 


y 7 


Mach line PA 








iat Plan form (ry-plane}. 
(b) Section уму, (rz-plane). 


FIGVRE 1.—8ketch illustrating coordinate system and the two degrees of freedom а and А. 
212483—54—56 
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In order to satisfy the boundary-value problem for фа, a 
convenient procedure is to start with the expanded form of 
the potential of a uniform distribution of doublets. Then, 
for a given power of the frequency of oscillation this poten- 
tial, as will be shown in the following analysis, can be 
modified so that, when integration is made over the арргорг!- 
ate region, the results satisfy the differential equation (1) to 
the given power of the frequency and satisfy the condition 
of tangential flow exactly. The type of doublet required 
is that with its axis normal to the plane of the wing. The 
potential of such a doublet may be obtained from the poten- 
tial of а source, located ш the plane of the wing, by partial 
differentiation with respect to the direction normal to the 
plane of the wing. Similarly, the potential of a distribution 
of the required type of doublets can be obtained from a dis- 
tribution of sources located in the plane of the wing. 

The potential at (х,у,2) due to sources located at points 
(2, 7,0) in region r (illustrated in fig. 2) of Ше ry-plane which 
satisfies the differential equation (1) may be written as 











ара o 
where 
- Ме Мо 
ധ്‌ cB? УВ? 
апа 
R— Cree 





FIGURE 2.—Sketch illustrating region of integration for the velocity potential. 
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Expanding the integrand of equation (6) into a power series in ш, collecting terms with respect to £, and differentiating the 
resulting integral with respect to z gives the expanded form of the potential of а uniform distribution of doublets, namely 


UE CA 9 [| | ug te. . „еи... (ада oe Та ћи... 1 ന 


т O2 
(ара... а...) аве (7) 
where 
a, m Oe (8) 
and 
(= (А (% а I A 
=m  "" (9) 
For convenience in the succeeding discussion and analysis equation (7) may be written in the following form: 
Ра д (—1)y" x) [| т-3 | 
v= 2а | [Рава + у >ы) аа (| зава (10) 


An interesting and significant property of equation (0: із that the coefficient of each power of £ satisfies the differential 
ғалы; 
equation (1) and has the form of а source potential with strength proportional to (=) е“. This property may be 


shown by writing the coefficient of 2% as follows: 


Any ( 1)%- =: а COME © 
Exe — 2) (м R) Tp vm +78) 


-(8y " [ е-ї#®-Ё® cos (% е) | (11) 





R 


А more general solution to equation (1) may thus be obtained by introducing properly chosen weight or distribution func- 
tions (denoted by £D,(£,9)) into the coefficients of Е" in equation (10). Let this solution be denoted by фі; then it 
can be written as 


а [പേ [Den മു മേ Б AG anf [рит нв ван | 12) 


m=? 


Examination of equation (12) shows that, at the surface 2—0, the potential ф is determined by the first integral expression 
but that both integral expressions may give rise to normal velocity. In succeeding steps in this analysis it is shown that the 
distribution functions .(Ет) in equation (12) can be determined so that the first integral expression taken alone will exactly 


Офа 


5 x) = “Уа, Also, any additional normal velocity that 
з= 0 





satisfy the boundary condition of tangential flow for фа; that is ( 


arises from the second integral expression can be canceled, to the required order, by consideration of additional doublet solu- 
tions to equation (1). Тһе problem of satisfying the boundary-value problem for the velocity potential фа may thus ђе re- 
duced to that of determining the appropriate distribution functions and additional solutions to equation (1). 

In order to show that the first integral expression in equation (12) can be made to satisfy Ше boundary condition for Фа; 
the coefficient 17" appearing in this equation is first considered. If the analytical expression for the coefficient 441 
(equation (8)) is multiplied by 2" and summed with respect to n, the result is identically Ve. This result may be shown 
as follows: 


Ра >) ағама Ра >) veer eter — Vagfszo 719 Уа (18) 


It is significant that this identity holds if only terms іп 9 to any given power are considered. For example, retaining 
only terms — ш to the fifth power gives 








и FECE gig e) (2 таа _ tOr? _ ve ERE) (-% кз, wirt iwr? 
Ye | (i-is “6 * 34 120) rot Fe (കി 4 12 J* 
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Next consider the normal velocity at 2—0 associated with ¢,, namely 


Examination of this equation and equation (13) reveals that the first integral expression on the right side of equation (15) 
yields Га exactly, provided that the distribution functions are determined so that the following integral equation is satisfied: | 


Ve 


M ПА 
Е Іш 25) | Dn) gd dq—Vaz (16) 





The kernel of this integral equation has the form of a steady-state doublet potential. The problem of determining the 
distribution functions for this case is therefore analogous to determining distribution functions for certain steady-state problems. 
The distribution functions for steady-state problems, at least for those involving conical flow, can be determined by a straight- 
forward process. The main details of this process are given in the appendix, where a method of solving equation (16) for a 
triangular wing is derived and the distribution functions required to derive the velocity potential for this wing to the fifth 
power of @ are given. It is to be noted that the method derived for solving equation (16) for a triangular wing may be 
generalized to apply to various plan forms and to problems of satisfying the boundary conditions for various velocity 
distributions. 

From this point on, the analysis is restricted to the derivation of terms of the expanded potential involving © to the fifth 
power. The method of deriving these first few terms is quite general and can be used to obtain as many additional terms of 
the expanded potential as may be desired. As previously pointed out, however, terms of the potential after the first few 
become very unwieldy. 

If the appropriate distribution functions are known for terms involving മ്‌ to the fifth power, equation (15) may be writ- 


ten as follows: 
Y= Vatw, : (17) 
where 














а. oO 2: 2200 gu ds = ees 1 4-4 Tz. 
к= — lim ааа at | |(1- i + മം intar EDR- മ മേ PDR- 


~ 


^ m 1 734 yv 
(поли тада Јер ED | (18) 


is the additional vertical velocity arising from the second integral expression in equation (12) involving @ to the fifth 
power. In order to maintain the boundary condition for фа, this additional velocity ഗം must be canceled. As pre- 
viously pointed out this canceling, to a required order, can be achieved by considering other doublet solutions to equa- 


tion (1). For this particular case consider relations ¢, and Ф. similar to à; (equation (12)) having the following forms: 


€ y д "m. a TE реса со — 1) -1 ny 21086 —2 » "€ 
втората te, | Рио рана + Бо (је) anf (Dat mem cas] (19) 


н-0 т-2 


== [а | Таз алад суа (С) а തടു | 195) 





The vertical velocity distributions 10; and иа, to the fifth power of 2, arising from these expressions are 


a ІМ | а: rris 8-4 m 227 “== P. 1053 ~ | NUT Pet 42 — 1 ra | or.) — 1 
ക lim EC (эн) f] | (iia T^ яр gt(io rau ts ) BD, :-(5- 5) BD, 5- 











Qe 2—0 Oz 

ധീ 1 c — AT га? 

в D p ayp во yy DR | ates к. 
— Ya. O /mW A „а 
бура а sa (17) Па із), то £D, А d£ d (20b) 


In [these equations the distribution functions Ор, Dy, and so forth сап be determined, as discussed subsequently, by the 
method given in the appendix so that 201 4-w, is identically equal in value but opposite in sign to и». When these functions 
are determined, the boundary condition for фа is satisfied by w,+%,1%.=Ve, which implies that the potential фа to the 
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fifth power of w is given by the sum 


фа=ф+ф + (21) 
to this power of wœ. 


Expressions for D, (n—0, 1, 2, 3) end D. (n=0, 1) are given together with expressions for other distribution 
functions D, (n=0, 1, 2, 3, 4, 5) in the appendix. In regard to the determination of D, and D,, it may appear necessary, in 
order to formulate integral equations for these functions, to perform the generally unwieldy integrations of the type 


lim Уз | [Бр mat dy (22) 
zo 02" Jr 


appearing in equation (15). In general, however, the information necessary for the determination of the functions D, and D, 
can be obtained, as is done in the derivation of the functions D, in the appendix, by examining the values of these integrals 


and their derivatives with regard to the parameter 8-4 &t some particular value of 6. 


Returning to equation (21) and introducing into this equation the expression for ൭ (equation (12)) and the expressions 
for ф, and ¢, (equations (19)), each to the fifth power of a, gives for the potential фа to this power of @ the following result: 


вота lim = | | јер, | (1—«а кы 0 YF R -(sxa- ശ്‌ ian) Gu — А) |+ 
ЕР | (ae LE Е Ее NE ET п) Rt san |- 3 dD, | (РЕНЕ 
(raw) [+ Pv ie ер (д “altel ады 

zia D (im 25 E Р (ജ്‌ 24% в edat D| (a+ zip |- 


sat | (2—2) а |да (а) Рада авер нне [60 |} а (23) 


Since it тау be shown, as in reference 8, for example, that 


Жара x | [^E dt Уат (г, y) 


T i0 Oz 























equation (23) reduces, at 2=0, to 


ത്ത "p 3 
€ =Val 14522 2% MR 2 ТАСЫ ти. 052 Ж) Deco Сс PEE = a Pe s (5-8 tote _ 





2 6 24 120 2 
БА ска) Пат > па ат. 57 24 pHa 9 i) Рух y EE 2 D (т, gt 6” E ie 
SS) Due nt (it its D) D (e 8) D арамы 
Zap ea) De egg Ри | (24) 


which, after the expressions for the distribution functions given in the appendix are substituted and the terms are 
regrouped, may be written in the following simple form: 


= Va O*—3? [Ao intei — во o ,B5y*) то а 13 + с В? yz) Huwes св? +018405) — о са 
ci 2? овча) (25) 
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The quantities A, с, and so forth in this equation are part 
of a group of quantities А, and പിപി which are 
defined in the appendix; these quantities are suitable for 
writing the potential and expressions subsequently derived 
for forces and moments in simple form. The quantities yj, 
Р, Q, and В; are functions only of the product BC (ratio of 
tangent of the half apex angle of the triangle to the Mach 
angle) and are tabulated for some values of 8C in table I 
and are shown plotted in figures 3 and 4, respectively. The 
quantities പം which are functions of Mach number M and 
the product 8C, may be evaluated for particular values of 
M from the values of Р, Q;, апа Е, in table I or from the 
plots in figure 4. 

The quantities A, and 4, are the same, as should be 
expected, as the parameters associated, respectively, with 
constant angle of attack and constant pitching of triangular 
wings of references 9 and 10. 

Expressions for the potentials фа and ф сап be obtained 
bv the method discussed for obtaining фа, or they ean be 
obtained to the fifth order of & by comparison and synthesis 
from equation (23). After simplification these expressions 
are 


фаст a Cry Ar — isc — 21383) — 401° — osB^ y x) + 
ге бег tH onp yz? 4- ci, 8595) 3-9 (ах + слову“ + 


сив —; ви (26) 
bah Ctx? у lAo ivre: — о (oar! + о у) + 
га? а Въ) Бебета сива сви (27) 
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At BC—1 or C= which is the condition at which the Mach 


lines from the apex of the triangle coincide with the leading 
edges of the triangle, equations (23), (26), and (27) reduce, 
respectively, to 


7 oan —2 
(де ау — DU (5 +7107 — 


ve: | 

(5—2M?)8'y"| + 7 3-8 M?)2— 
(1 —2M?yg?y*z] Fran (833/*--414342--63)*-- 
42(23/3 — 9M? — 3)8'y221--(8315—36344-63)85] — 

10) 73 л. Ру н „= 
31558 ың ТЗ M4638 3--297):4-0(46M* 
253 M? —297)8*5z:-E3(8M*—44 M? +99)pty*z]} 
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297)254-2(118.M*—1221M? —297)8*y'z* — (30) 
| Equations (28), (29), and (30) can also beobtained by going 

(168M*—396 M?—297)8*y'z] d ЗЕ (фикс: (29) | to the limit $C—1 in the expanded potential for triangular 
wings with supersonic leading edges (see, for example, refer- 

ence 2 for results to the third power of 9) or by integrating 
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the expanded potential of & unit source over the region oc- 
cupied by the wing with 8C—1. Equations (28), (29), апа 
(30) therefore serve as a check on the results given in equa- 
tions (25), (26), and (27), respectively. 

For values of the product ВС such that 8'C*«1, equa- 
tions (25), (26), апа (27) reduce, respectively, to 


dar V a, СЗх:— у: (81) 
фа, C* rt — y! (x — хо) (32) 
фа==ћђ, Сх? у? (33) 


These expressions are the counterpart of the potential associ- 
ated with triangular wings of vanishingly small aspect ratio 
in steady flow (reference 11) and may thus be regarded as 
the potential associated with a small-aspect-ratio triangle 
oscillating in pitch and vertical translation in either sub- 
sonic or supersonic flow, provided that s remains finite. The 
condition 8C«&1 is obviously satisfied for all values of C 
when 1f=1. However, as М->1, the value of w becomes 
infinite and the expanded. potential, as treated herein, be- 
comes Meaningless. 
FORCES AND MOMENTS 

As pointed out in the introduction, the velocity potential 
for the triangular wing can be used to calculate the aero- 
dynamic forces and moments for other plan forms that can 
be formed from the triangular wing by cutting the trailing 
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edges so that they lie ahead of the Mach cones emanating 
from their foremost points. Sketches of different plan forms 
thus obtained are shown in figure 5. 

The force апа moment coefficients desirable for most flutter 
calculations are those that yield the spanwise variation in 
these quantities or chordwise force and moment coefficients. 
These coefficients are obtained by integrating the pressure 
difference along any chord for the forces and the pressure 
difference multiplied by а moment arm for the moments. А 
convenient procedure in deriving these quantities is to intro- 


duce the reduced frequency parameter bek and to employ 


the variables 2, y, and хо in а new sense as nondimensional 
quantities obtained by dividing the old variables by the 
maximum chord 2b of the wing. 

'The pressure difference between the upper and lower sur- 
faces of the wing 18 
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In equation (37), for example, the quantity (L,--4L4) is the 
lift-force coefficient associated with vertical motion of the 
wing. The real part Г. is in phase with the vertical position 
of the wing and the complex part L, is 90? out of phase with 
this position. Similar definitions apply to the lift coefficient 
Са КЕ) associated with pitching motion and to the moment 
coefficients (M, +iM:) and (Мұ--4М). The imaginary or 
out-of-phase terms determine the aerodynamic damping 
associated with different wing motions. 

Although the expressions for the components of lift-force 
апа moment coefficients in equations (38) and (40), respec- 
tively, are lengthy, they may be quite easily evaluated with 
the aid of table I or the graphs in figures 3 and 4. 





ടു 
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(8) Plan form А. (c) Plan form C. 
(b) Plan form B. (d) Plan form D. 


FIGURE 5.—£ketches (ustrating different plan forms for which the force equations (38) and 
moment equations (10) apply. 


DISCUSSION 


SAMPLE CALCULATIONS 


In order to give some indication as to the general nature of 
the spanwise distribution of the different components of lift 
and moment coefficients, equations (38) and (40) have been 
evaluated at different spanwise positions y for plan form A 
(fig. 5(а)) for the following set of conditions: BC=0.5, 
С=1.0, 2--0.6, M—45/4, and k=0.1. These sample re- 
sults are plotted as functions of spanwise position in figure 6. 
The spanwise variations of the different components of lift 
force are shown in figure 6(a) and the corresponding varia- 
tions of moment coefficients, in figure 6(5). 

In figure 6 note that, for the particular set of conditions 
for this example, the maximum numerical values of the com- 
ponents of moment coefficients Mi, Ma, and Af; are positive 
and are near the tips of the wing, whereas the maximum 
numerical value of the component Af, is negative but is also 
near the tips. It may also be noted that the integrated 
(in spanwise direction) values of the components of moment 
coefficient, or components of total moment coefficient, would 
in each case have the same sign as the maximum value of the 
corresponding component of section moment coefficient. 
This result is not necessarily true in general, because changing 
some of the parameters involved in the evaluation of the span- 
wise distribution of some components of both foree and 
moment coefficients may change the distributions signifi- 
cantly from those shown in figure 6. | 

Тће fact that the total component of moment coefficient 
M, is negative in the example just discussed shows that, for 
the conditions of the example, this term would not contribute 
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їз) Lift force. 


Fi uEX 6.—Spanwise distribution of components of lift force and moment coetfictents for 
14+0.8. 8С-0.5; Cm1.0; k=0.1; Ме VE | 


to the aerodynamic damping but, on the contrary, would act 
as a source of energy for the oscillating system. This cir- 
cumstance is significant since it leads to the possibility of the 
single-degree-of-freedom torsional instability discussed in the 
following section. 


UNDAMPED TORSIONAL OSCILLATIONS 


The wing plan forms discussed herein, like two-dimensional 
and rectangular wings, exhibit the possibility of undamped 
torsional oscillations for certain ranges of Mach number М 
and location of axis of rotation г. This fact is borne out, 
as indicated in the preceding paregraph, by considering the 
integrated (spanwise) value of the component of damping 
moment Af, associated with pitching or torsional motions. 


4 
Ma 
0 
“ы 
(b) Moment. 
ЕЕЕ 6.—Concluded. 


The main results of this phenomenon can be obtained by 
considering very slow oscillations so that only terms in 
equation (40d) for M, involving the reduced frequency Е 
to the order 1/k need be retained. Іп this case, 
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For the triangular plan form (fig. 5 (a)) the integrated value 
от total component of moment coefficient is 


17,701 LM) ASESQM?— 1) А,+ 


28*k 
Aroldo (2 4—1) A] t 42 CM? —1)4,) (42) 





In general, the condition of torsional stability or instability 
depends on the sign and magnitude of Af, in equation (42). 
Positive values of 27, indicate stable conditions and negative 
values indicate the possibility of torsional instability. 
Between the stable and unstable conditions—that is, when 
AM, vanishes—a borderline state of unstable equilibrium 
separating damped and undamped torsional oscillations exists. 

The ranges of values of Mach number М and location of 
axis of rotation го for which Af, vanishes for some selected 
values of C=tan є are shown plotted in figure 7. The 
regions inside the curve in this figure indicate instability. 
The dashed curve, on which some of the solid curves ter- 
minate, represents the locus of values of Af and 2, for which 


M, vanishes for the whole class of triangular wings with 
‘supersonic leading edges; that is, for triangular wings where 


ВС>1 (see fig. 5 of reference 2). It will be noted that JM 


vanishes for values of 2, ahead of the root chord position. 


It wil also be noted that, as the vertex angle e—tan^! С 
decreases to 30°, the range of values of Mach number for 
which M, vanishes decreases sharply. 

In conclusion, investigation of equation (42) shows that, 
for а given value of the reduced frequency Е, Mach number 
М, and location of the axis of rotation zy the magnitude 
of the damping coefficient 24, generally decreases as С de- 
creases and, consequently, torsional instability is less likely 
to occur with slender triangles than with wider triangles. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
LANGLEY FIELD, Vå., June 19, 1951. 
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FIGURE 7.—Curves showing ranges of Mach number М and axis of rotation ze for which the aerod ynamio torsional damping moment vanishes for some selected values of angh e. 


APPENDIX 


METHOD FOR DETEEMINING DISTRIBUTION FUNCTIONS 


In this appendix a method of obtaining the distribution functions ts developed in detail. Expressions for the distribution 
functions required to drive the velocity potential to the fifth power of the frequency of oscillations and & set of functions 
useful in writing the expression for the potential in concise form is given. 

As indicated in references 9 and 10, where the triangular wing is treated for constant angle of attack and for constant 
rolling and pitehing motions, а convenient form of the velocity potential corresponding to а distribution of vertical velocity 
proportional to z" is 





ps tide 
Е, TT Al 
али Eus an (2 — E — By — c E — 822? = 


where є is the slope of a ray passing through the vertex of үн: wing (п=<Е), & із the least value of Е that causes the 
denominator R in the integrand to vanish, and F,(c)= =7 t DEn) is the distribution function that ја to be de- 


termined so that 





(വു) == (A2) 


In equation (А1), the integrand is noted to be singular at the limit &= ё. А form of the integrand which avoids this 
difficulty is obtained by making the following change of variables: 


q Би (1— 8726) 


N= D 19൦ (6- cy +5 G= 65൦) 








(43) 
т | 
т 
r=cosh™ 2—2 
where т is the new variable of integration. With these substitutions equation (Al) becomes 
10 Fata) $ сЕ 
Фа x 3z Caa, (q—N cosh T) 2dr (А4) 
апа the corresponding expression for ш„=( 90) › which is an integral equation for Е (с), is 
QVE Иги 
zx 
l.. "e Fa cosh Ар Ра atil 
C=- lim — amo |, А (cosh г) cosh т т + 
а?) 1 __Едајда _ ES (5-соз (അ cosh ) cosh d 
а— (1— gigs оз? |, 1 ү Т N T тат 
=>" (А5) 


In this equation the value of и, obtained by performing the indicated integrations and then going to the limit 2=0 is the same 

as the value that would be obtained by first going to the limit &nd then performing the integration but neglecting singularities, 

pointed out subsequently, that arise when the value of с approaches the value of 0. Making use of this fact reduces the 
579 
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calculations involved and simplifies the integral equation for Ка. Thus passing to the limit gives for equation (А5) 


Mg = BuU ce И ден? Aw UA Ма" (2 —cosh г)” cosh т dr 


L^ | (A0) 


where №, is the value of № аё 2—0. The requirement that the normal velocity be proportional to 2" and independent 
of y implies that all derivatives of equation (Аб) with respect to y or 6 must vanish or that Е (а) must be so determined 
that the final value of the integral in this equation be independent of 6. Тһе requirement that all derivatives of equation 
(АВ) with respect to 0 vanish leads to other equations for Е, and, after n+1 such differentiations, the equations 
acquire forms for which solutions are known. Тһе value of Ғ, can then be determined from these known solutions 





Е 
by evaluating each of the derivatives 2 T at any arbitrary value of 8 in the range —C <8 <C as follows: 
The Ath derivative of w, with respect to @ (equation (А6)) gives the following integral equation for Ру: 
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= Bo-+cosh т -X) 
аш, Мо cosh т-а No 


(1) от 
(No (а — М» cosh т)" 


ച (n 1112൦ (А8) 
X (– Мр“ (а — № cosh т): 


Thus the integral equation corresponding to the n4-1 derivative of w, is 
гуё 








മട (o дей [ട്ടം 
@@ її т -е NP — 820595 0 
Иља T АЕ А ш. | 
== - _с(@— с)” да —0 (А9) 


Further differentiation of ഡം leads to other equations involving integrals similar to that in equation (А9) and, as will be seen 
subsequently, is not necessary for the determination of Ра. The singularity at с--0 in equation UM ls & result of going to 
the limit z=0 after equation (А5) and, as previously implied, is to be ignored. 

Consider the following equation similar to equation (A9): 


( С Дао) е o (А10) 


-c(8—c) . 


It is known by analogy with problems in incompressible flow and may be shown by direct substitution and reduction that 
this equation is satisfied for any value of n (п=0, 1, 2, . . .) by the function 
(о) =ү02— 0? (A11) 


Е, 
. This function also satisfies equation (А9) for any value of n and satisfies equation (AT) for all derivatives T E: beyond 
the nth derivative. А more general form of solution to equation (А10) may be shown to be 


f= oOo? | (412) 
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where m is an integer and m €a. The validity of this solution follows, since 


[കിട 33 a 0-2 (A13) 
and 
o".CiI—. в. т!6е- Елан Д 


(8— а = (т —г) 171 (0—0)" "+8 (Ам) 
Each term of the summation in equation (А14) 15 found to have the form of the integrand in equation (А10), from which 
the function (A12) is concluded to satisfy equation (А10) and equation (А7) for all values of K>n. 

From the foregoing discussion and consideration the distribution function Е (а) may be uniquely determined in terms 
of expression (А12) as follows: 

Consider the expression . 


Едо) --Вл+-Вл+ cee TB. f. (A15) 


where the coefficients By, Ві, . . . B, are constants that are to be determined. 

Each term on the right of equation. (А15) is noted to satisfy equation (АТ) for all derivatives beyond the nth and a 
total of n+1 parameters is to be determined. If expression (А15) в introduced into equations (А6) and (A7) and the indi- 
cated integrations in equation (Аб) and in equation (А7) for K=1, 2, ... are performed, n+1 linear equations are 
obtained in By, В, . . . B, from which these constants шау be determined. . 

The integrations with respect to с in equations (AG) and (А7) are in general difficult and tedious to perform; 
however, as previously pointed out it is only necessary to perform the integrations for some particular value of 6 in the 
range —C«6« C. The integrals have their simplest form when 6—0 and the integrations can be made for this value of 0 
by reductions and use of formulas in reference 12. 

The functions D, and Б, тау also be determined by the method discussed for determining the functions D,. The 


conditions to be satisfied by D, and D., however, are determined by the value of w, equation (18), and its derivatives 
with regard to 0 at 0—0. B 

The distribution functions D, (n—0, 1, 2, 3, 4, 5), D, (n—0, 1, 2, 3), and D, (n—0, 1) calculated by the foregoing 
method are as follows: 






































_ 4 Te | TAT 7 ту AB 
275൧൦0 Е--1 0,-(2,+2, g «4C Е-т 
Рът А „СВ D, (4, 8 кр 022—1? 
2 2 да а 
р (Art A 52.) CEP D,-( 2+4, = БАН ма, Cg 
2 — M — _ Rm. _ д4 
Әу-(.4,--4, BL) „Сет 0,-(2,+2, 2 +4, ET) уС?2— 1? 
Dy=( Act Ar GP p” т P AP B x) сет D- (2-4. e 4-4, E a fr), Ce 
22 4 „4 NOS == = = 
=( A+ Au Рада e \ ORF D;— (4, ന A +A, E ET Cg 


where, with the notations 


| 2 PHO m coste S. n+l 
റ 722)8 | Oo de lim 25 д? (oer лга (-сөзі r) cosh т dr 
0 





тх" J-e (0— Bh? 967 Jo 





= 
== 
zm 


nci 2 +C gm iC 22 ОР | созћ-:— 
Tp а | =т= de lim 557 35 х (а —N, cosh ന [2g —(n.4- 4) №, cosh 7] sinh? т dr 


and 


д? 


26? =>“ NXg — №, cosh ന)": [4g — (n-]- 6) №, cosh т] sinh* т dr 


= 2 
т.а |. с m. gi g^ 41— 1— 880° do lim 
T -с 


тв 


( Note that the above expressions will vanish identically unless p and m are either both odd or both even.) 
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Ше A’s are given by the following equations: 





1 
A= - 
0 We, 
1 
А у= 
1,0 


A Po +A P= 
യം 
АБА, 1 

x cd 


ВИ AW: А.И. — 
А, o AM reat 
As Wiot AW А Wt «-- 0 


АЙ РАИ РА 1 
АИ otA V3, 1-4,1 W= 
Wiot А о Ws + А Wi 0 


AoW? + 4, W2 „== AW. 
A, 


УЗ tA, Wi. mes 

AW? + А 9 =A, Wo 

A, W? + А W = AW? 

А.И t А Ра + Ае = А ТУ о + А Ws 


A3 d AW + AW А, Ио А.И, 
А-А ва AW AW, +A Wis 


_ (27р --81р + 12 po°)F" + 


A, Wot Ав БАИ АЙ РАЗ, 
A те 2, +A W = A, EA, 
: AST: а AW, «= А, H^ о ASI. t 


МА 


1 


= - EN ej? e 
| 
| 


AW, NNI m4 


AWe o АЛО ҚАНҒА = — 
AW АН, А, — 
നി БАЛ, БА = 


where, with 8C = po, 


"eor E° 
ра PEE rü-—20)E'. 
" (1— po’) 
gro (Spd —3 ps) P" (2— 10 po? број E" 
2,0 | Hi — 55 
wa -2р (മേ оо) Р” 
2.3 8105) 
а 2—26. W, 
we _ (бре —Зројћ  +(6—102»2 + 220 Е” __ 
M 1— ро)? 


(6 — pan i ha 


Wo _ (56р —92р*-Е72ре° —20во)Ё' + +(8— — 117 po’ +202 pot— 149 p +40 po) E^ 





8(1— po. я 


W$ 
6(1— po)” 
Ww? рио) (Вр ET oi Bee)" | 
i 68*(1— ро) 
2 
Wis ЕК 
Wi, നനനനത്നനത്ത്ത്ന്‌ 
(1— род? 
88*(1— ро) 
We EPET EPAF (ра 142൦-20. — 
881 — po) 
2 
И = т Wia 


(609 —30 p —30 po? + 6 pF’ —(12 —39 po! + 6 o9 — 39 p+ 12004" 


Wi മപ 


2(1— po * 


1 70 А, 70 = ന +4, 2, 


а БАЛОНА, 


^ + АНА НА," 
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үз же Әрі — во 51 ൧൦7 — (бр —242 0! t 54 po°+ 12 р) E" 


28%1— pc Р 

де: 2 

ae Wea 
Wi se “Wi, 

422 (16808-97 
И + = бёр 27 6 po° + 21600 — 60р) + (24൦351൧൦ 4- 606 p, 447 p 2.1. 120)Е” 

(1—20) | 

ә __ (4002 —833 po +994p—553p + 12009) Е" --(40 —859 ре + 19102, — 211 5ро 11362, — 240 p01) E 

Ke о 0 ро +994 ра 558р 4- 120009) Е +(20— 859 0124-19106, —2115p454-1136p,5 240p 10) E” 
0 


10(1— рој“ 


6 ғ 
ит __ (5500+ 7400 —ро YF” —(5p 108р +17 ро —2po9 E" 








408*(1— po”) 
> 208? 
їр = m We, 
ലം 11 2ро 49 ро бреје "7 (12—4Тр +—22 pot—127 ро 
io 0 po 2 Гро 72р — 165p) М” 
Wi _ (8 ро“ —2р04-1310; —4pj4 9) F" —(6 ре — 8120 
ро'+108р2 0 +53 ро — 8 po Е“ 
то сети. Ро ро ЗЕ 
anes 
ПЁ ру 
" 68? 
Wie Wi. 
Wi _ (24 po? + 609 p,* — 462 pa +273 ра — 60 p.*)F’+(120— 
—561р  -+102 — j 
ПОР? po 102 6 po* 969 204-48 p," —48 pol E 
T ഗ്‌ 
| 2(1— po) 
Wz ഗ്‌ 
(1— po) 
വ .=— (3 poi рођ! —(4ру —2 po) E" 
6(1— ре)? 
T കി 
3 ൧൧൦൧7൦: 3000) —(19 ро — 17 pot + 6 pE” 
24(1— po 
pm (3 po + броје —(7 po + po) E" 
248'(1— — ро)? 
13 Вог ps) F" —(2+7ра — po ] 
2(1— psy 
ന്‌ (pot tT po) Е" (pd! — 7 po —2 po) E^ 
2(1— ре)? 


Wi „— — 8 Бро Зао) UETA] 
(1— ра) 
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120(1 — po 


To _ (15 pof 84 po — ро) —(38 ppt +12 po —2 pF) E" 
1208'*(1— po^? 








Tre 8 КЗЗро 3 14p0 --роЈе'—(6+46р о —бре "гр 9E] — 00 ha _ ПЦ У 
3.0 6(1— p 

Te (009-7440) —2р0)" (30൧൦-34 2100-40). 
= 6(1 — pot 





rs g=— } За 

= 2 

У ~ ке 

Tr 8B lee ро) — гр Е| | | u "NR s 
i 2(1— po 

— 4 

Toe ИЗ: 

= 2 

Tm = We, 

== ыы 

ന്നു 


Combinations of these functions useful in writing the potential in concise form are as follows: 


R 
= Pee Se 


where 

Ре Ap—~ d; @,=0 К,=() 

P,= А,— Л» Ооғ- 0 250 

(൨൭. Ao 

= (Ap—2A,+ Ay) t= В.--0 

Pa Q-—5 S R0 

Ре (Ao Bet 4-8.45— 44) 0---5 (Asi) Ry=0 
1 1,— = 

2൦ (84,— 4) Q——5 (A1— Aa) R,=0 

P=} (A,—2 4:4 A) y=? IE R,—0 
1 A. 

Рұ--2 (24,--А) wz R,=0 


E 


боо 


4. 


5. 


МАСА. Вер. 872, 1947. 
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Py (4--441--6.А,-44А,4-.4)) 





Pu ry Ы 
| 
P» 120 (A p— 5.3; 3-10.45 — 10.4, 4-5. 4,4 — Ay) 
1 
Роб (104;— 10As+5A;— Ain) 
(ИЕ 
14-7 190 (5.A5 — Ал) 


1 
Pus (41—3.3, -3 ചി. 44) 


1 
Р 15 6 (3 4 —3 Aa БЕ Аз) 


ds 
P= 6 


Pa (A5 —4 Ae 6.4, —4 А,+ Ag) 


1 
Р "=g (Ар—44; + 6.4,- s) 


Pa (4:4 —4 Ag) 


Ф--:(2-22.-2) 


eS ur 
Qi——1 (A; —2:4,-- A.) 


ഹ 

== -1% (4-34,-4-34,-.4;) 
P vee — Ви. 2 

Qu——15 (4,—34,+34;—2%) 
1 Зе ә 

Qu——15 (3.4,—A,) 
im вез 

@== 2 (А.-А) 


Qis= -i (45-542 


ды--1(2.-22,4-2) 


Q= -i (A,— 24, As) 
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